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A self-consistent set of equations for the one-electron self-energy in the ladder approximation 

r^ ' is derived for the attractive Hubbard model in the superconducting state. The equations provide 

an extension of a T-matrix formalism recently used to study the effect of electron correlations on 

normal-state properties. An approximation to the set of equations is solved numerically in the 

0^ ' intermediate coupling regime, and the one-particle spectral functions are found to have four peaks. 

This feature is traced back to a peak in the self-energy, which is related to the formation of real-space 

bound states. For comparison we extend the moment approach to the superconducting state and 

f-H I discuss the crossover from the weak (BCS) to the intermediate coupling regime from the perspective 

Q ■ of single-particle spectral densities. 
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I. INTRODUCTION 

High-temperature superconductors display a wide range of behavior atypical of standard band-theory metals. In the 
superconducting state these materials become extreme type II superconductors, with a very short coherence volume, 
which one might take as an indication of tightly bound pairs. Particularly interesting in this respect are recent 
experiments showing a pseudo-gap structure in the normal-state density-of-states of underdoped Bi2Sr2CaCu2 08+5 
that persists almost up to room temperature |||,g|. These features suggest that correlation effects strongly affect the 
properties of such materials. One of the simplest models featuring superconductivity and allowing a systematic study 
of the effect of electron correlations is the attractive Hubbard model. Although this model is unlikely to provide a 
microscopic description of high-temperature superconductivity, it is likely that it reveals the effect of correlations on 
measurable properties. 

In a previous work, the effect of electron correlations on some normal-state properties of the attractive Hubbard 
model was studied using a self-consistent T-matrix formalism, going beyond simple mean- field treatments M. It was 
found that for intermediate coupling strengths the attractive interaction gives rise to large momentum real-space 
bound states with energies below the two-particle continuum and with a pronounced effect on the spectral properties, 
namely a splitting of the free band into two, one of them associated with virtual bound states. Also, a bending in the 
static spin-susceptibility was observed for temperatures just above the phase transition. 

Thus, it is of significant interest to study the effect of these bound states also in the superconducting state. To this 
end we derive, in Appendix A, a set of equations for the electron self-energies in the ladder approximation including 
anomalous diagrams. This is done in a general fashion using the technique of functional derivatives, well-known from 
the textbook by Kadanoff and Baym H. To facilitate a numerical analysis, the anomalous part of the equations is 
approximated in a somewhat simpler form, which is a reasonable approximation for the parameters of interest. These 
equations reduce to the familiar T-matrix equation in the normal state and provide an obvious extension of the work 
embarked on in Ref. Q. In Sec. || we present a numerical solution using the fast Fourier transform (FFT) technique 
and discuss the physical meaning of the results. 

To gain further insight into the electronic properties of the superconducting state we compare the results with a 
generalization of the moment approach, allowing a description of the diagonal and off-diagonal single-particle spectral 
function in the superconducting state, concentrating on the weak and intermediate coupling regimes. This formalism 



is developed in Appendix B, and a comparison is presented in Sec. [II. 

Our main results include (i) the appearance of four excitation branches in the one-particle spectral function, with 
symmetric dispersions around the chemical potential, and (ii) that the low-energy behavior with respect to the chemical 
potential is in qualitative agreement with BCS, whereas the high-energy features are not accounted for in BCS. 

II. T-MATRIX EQUATIONS BELOW Tc 

The Hubbard Hamiltonian is defined as 

H = -t ^ cl^ci-a + U^ni^nii- ^i^ni^ , (I) 

<ll'>a I la 

where cj^icia) are creation (annihilation) electron operators with spin a. The number operator is njo- = cl^cia, t a 
hopping matrix element between nearest neighbors I and I' , U the on-site interaction, and /i the chemical potential 
in the grand canonical ensemble. Here we consider an attractive interaction, U < 0. For a review of the attractive 
Hubbard model, see Micnas et al. [^. 

In Appendix A we used functional derivatives to derive self-consistent equations for the self-energy in the ladder 
approximation in the superconducting state. Here we present a numerical solution of a simplified form of the derived 
equations, where the superconducting order parameter is only taken into account to lowest order — appropriate when 
it is significantly less than the bandwidth. The use of a ladder approximation means that only two-particle scattering 
events are described and restricts the validity of the scheme to dilute systems. Thus, we approximated the self-energy 

byffl 

Tr^ -r'^ - /^ '^22(a;,a;')T(x',2;) A{x)Six - x') \ ,. 



A*(a:)(5(a;-x') Gii{x,x')T[x' ,x) 
where x = {I^t) denotes a space-time coordinate and the T-matrix in Fourier-Matsubara space is given by 

T{q,iem)^- T^-, — ■ — 7 (3) 
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x(g,i£m) = ^^7^ ^ G22(A; - g, iw„ -i£„i)Gii(fc,zw„) , (4) 



where iSm and itOn are bosonic and fermionic Matsubara frequencies, respectively. Another approximation approach, 
taking the order parameter into account in a higher order than here, can be found in Ref . Q . 

In this formahsm the two-particle Green's function, which is also the local-pair propagator, is given by H 

G2(k,z£„) = . (5) 



Green's functions are determined self-consistently using Dyson's equation, Eq. (A2), and the order parameter is 
determined by A (a;) = J7Gi2(a;^, x). Technical aspects of the numerical solutions using the FFT technique have been 
described in detail in Ref. [B|. 

The numerical simulations were performed in two dimensions (2D) for U/t — —6, a density of p = 0.1 and a 
temperature of T/t = 0.05, which is far below the transition temperature. Here one must note that, owing to the 
Mermin- Wagner theorem, no phase transition is expected in 2D in systems with a continuous symmetry, and the 
formalism employed is too simple to describe a Kosterlitz-Thouless phase transition. Nonetheless, it is observed that 
the formalism does yield a phase transition, in agreement with the results above the transition temperature, where a 
divergence of the T-matrix was observed with decreasing temperature. 

In Fig. n^, we present the diagonal one-particle spectral function, ^(k, w), defined by 

ACk,uj)= lim -~-Im\Gn(\<i,uj + iS)], (6) 

We see that four peaks exist for every momentum along the diagonal of the Brillouin zone, a situation we will 
examine in a moment. The valu e of th e order parameter is A/i = 1.04, which is within the region of validity 
of the approximation, Eqs. ( |A13| )-( |A17| ). The transition temperature can be found as the temperature where A 
vanishes and is T^/t — 0.2, compared to the BCS transition temperature T^^^^/i = 0.66. The zero-temperature order 
parameter, critical temperature ratio is now 2A(T = 0)/A;bT'c = 10.4, quite different from the canonical BCS result 
of 2A(T = 0)//cBrc = 3.52. 

In Fig. H, we present the off-diagonal one-particle spectral function, i3(k, w), which is defined as 

B(k,uj)= lim Im\Gi2(k,uj + iS)] . (7) 

We also observe four peaks in the off-diagonal one-particle spectral function which arc at the same positions as the 
ones of the diagonal one-particle spectral function. In Sec . Ill we will see that symmetric poles around the chemical 
potential are a consequence of Dyson's equation [Eq. ( |A3[ )]. 

To explain the occurrence of these four peaks in the one-particle spectral functions, let us study the real part of the 
self-energy. This is presented in Fig. 0. E xam ining the one-particle Green's function, G(k, iwn), coming from Dyson's 
equation, we find it is given by [see Eq. (|A3|)] 

^/, ^ ^ z-H£k + £(k, -z) 

^ ' ^^ (^ - £k - S(k, z)){z + £k + S(k, -z)) - |A(r)|2 ^ > 

using S22(k, z) = — Sii(k, —z). We wish to study the roots of the denominator of Eq. (||). For example, in the BCS 
case where the self-energies are zero, we obtain two solutions at u; — i-^/e^ -I- |A(T)|-^. For k = (tt, tt), we obtain 
£k/^ ~ 8.4 because e^ is measured with respect to the chemical potential, which in our case is —4.438 (see Fig. |^). 
Then, the two roots are located at ~ ±8.5. These two roots are going to remain for the interacting case, as we will 
see shortly. The positions of the poles are determined mainly by the real part of the self-energy, and the imaginary 
part contributes essentially only an additional lifetime effect. Thus, we assume for simplicity that the imaginary part 
is zero. We then have to find the roots of the equation 

(c^-ek-S'(k,c^))(c^ + ek + S'(k,-c^))-|A(r)|2 = , (9) 

where E'(k, w) is the real part of the self-energy. A rough estimate of the roots at k = (tt, tt) of Eq. (^) yields the 
values uj/t « ±2.6, ±8.5. We see that the two previous (i.e. BCS) solutions remain and that there are two other 
solutions symmetric around the chemical potential. These excitations correspond to electrons bound into local pairs. 



Figure ^ shows the imaginary part of the two-particle Green's function, — ImG2(q, ^)- We observe that for zero 
momentum we find the Cooper resonance, i.e., two symmetric peaks around the chemical potential. However, the 
Cooper resonance is different above and below Tc. Above the critical temperature the Cooper resonance is sharp 
and close to the chemical potential, whereas in the superconducting state the resonance is less sharp and the two 
peaks are separated by an energy of w 2A. For larger momenta, we find a peak that yields a narrow band (around 
uj « — 2/i + U). This narrow band is associated with pair states, which is a feature also observed in the T-matrix 
calculations above Tc B. The pair states give rise to a narrow band in the self-energy, which in turn gives rise to a 
band of bound electrons in the single-particle excitation spectrum. For this value of the interaction, the band crosses 
that of the BCS quasiparticles, causing a hybridization and correspondingly a split into four excitation branches. 

III. MOMENT APPROACH 

To further the understanding of the results, we have generalized the moment approach Q to the superconducting 
case. The formalism is derived in Appendix B, and here we present numerical results and compare them with the 
T-matrix results from above. 

For the T-matrix formalism, we presented numerical data for U/t — 6 above. This interaction strength was chosen 
because it provides a clear picture of the essential physics. Unfortunately, the moment approach has a limitation for 
higher coupling strength, as will be discussed below, and we found it necessary to restrict the coupling strength to 
U/t = 4 in the following. 

In Figs, g and|g we show the dispersion of the quasiparticles and their weights in the normal state for U/t = 4, p = 0.1 
in 2D. It can be seen that the attractive interaction splits the free quasiparticle band into two bands. Examining the 
weights, we identify the unpaired electrons as being the species of electrons with dominating weights — residing in 
the lower band for low momenta and in the upper band for large momenta. With increasing U the bands separate 
more, and for some critical U only paired electron reside in the lower band, whereas only unbound electron reside in 
the upper one |10|] . Retaining this set of parameters for the moment, we calculated the dispersion of the quasiparticles 
and their weights in the superconducting state usingthe moment approach outlined in Appendix B. 

The results for U/t = 4, p = 0.1 are shown in Figs. and g. We obtain four bands instead of two in the normal state, 
with a gap around the chemical potential. The lowest band has negligible weight (Fig. @) and is irrelevant for the 
diagonal spectral function, although it makes a finite contribution to the off-diagonal spectral function and thereby 
to the gap equation and Tc. An essential point is that the band wherein the chemical potential was located as well as 
the weight are divided into two, leaving the upper band practically unchanged. For this reason, a BCS treatment still 
describes the behavior around the chemical potential reasonably well, but misses the additional excitation branches. 
These affect the off-diagonal Green's function and hence the gap equation, Eq. ( |B36| ). Indeed, the BCS transition 
temperature T^'~^^/t — 0.37 is reduced to Tc/t = 0.19. Thus, whereas the BCS transition temperature increased 
without bounds with increasing coupling strength, the real transition temperature is reduced compared to that 
resulting from the increasing pair mass. Indeed, for infinite coupling strength all pairs would be localized and the 
transition temperature would go to zero. These features characterize the crossover from weak to intermediate coupling. 
In Figs. and g we have also included results from numerical T-matrix calculations for the same coupling strength 
and find quite a good overall fit. The discrepancy in the dispersion for large momenta is also found in the normal 
state and can be attributed to the neglect of the fc-dependence in the approximation for the band-correction term in 
Eq. (pTp. 

For increasing U we find that the chemical potential moves into the correlation gap and that the model becomes 
an insulator. The formalism used here yields Tc = 0, which is correct in the sense that we no longer expect a Fermi- 
surface instability for large U. Instead we expect to have a scenario of Bose-Einstein condensation of preformed 
pairs, yielding a finite transition temperature. The description of this phenomenon is beyond the current formulation, 
however. Considering that the moments describe the transition to the strong-coupling regime well in the normal state, 
we expect that this will also be so would we satisfy more off-diagonal moments. For instance in our approach all odd 
off-diagonal moments are zero. This, however, is not true for the exact moments. We thus believe that these nonzero 
moments must be satisfied in order to describe the crossover to Bose-Einstein condensation. 

IV. CONCLUSIONS 

Using a functional derivative technique, we have generalized the particle-particle ladder approximation to the 
superconducting state, providing an obvious extension of previous work in the normal state. The resulting gap 
equation is a complicated, nonlinear integral equation, whose solution is difficult to obtain. Recently, Bickers and 



White evaluated the T-matrix eigenvalues, which is equivalent to looking for normal-state instabilities that signal 
second-order phase transitions |11[| . In contrast, the equations derived here cover the entire range of temperature, 
both above and below T^. 

Another calculation was carried out by Haussmann [|l2| , who developed a perturbative formalism for the vertex 
function starting from diagrammatic techniques, including both particle-particle and particle-hole channels in the 
Bethe-Salpeter equation. He ultimately neglects the particle-hole channel, the vertex function is transformed to the 
scattering T-matrix, and he treats superconductivity in a mean-field fashion. Thus, no direct comparison with our 
equations is possible. 

A somewhat simplified version of the equations derived here has been solved numerically using the FFT technique. 
The main result is the appearance of four peaks in both the diagonal and off-diagonal one-particle spectral functions, 
which can be explained by a narrow band in the self-energy, related to the formation of real-space bound states below 
the two-particle continuum. 

We presented an extension of the moment approach to the superconducting state, providing an approximate ex- 
pression of the single-particle spectral density. Although lifetime effects are neglected and correlations effects in the 
band-correction term have been treated approximately, our approach provides new insight into the crossover from the 
weak to the intermediate coupling regimes. 

For the specific case of U/t = 4,p = 0.1, we have shown that the single-particle spectrum in the superconducting 
state exhibits four symmetric branches with respect to the chemical potential. The low-energy part (with respect to 
/i) resembles BCS behavior, whereas the high-energy physics is strongly modified. This led to a reduction of T^ in 
comparison with the expectation from BCS. 
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APPENDIX A: DERIVATION OF T-MATRIX EQUATIONS 

It is convenient to define the Nambu-Green function 



^^ ' ^ \ {Trcl{r)c\.,{r')) (T.ct,(r)c.;(r')) ] ^^'^ 



where T^- is the time ordering operator. 

The Dyson equation for the self-energies in Fourier-Matsubara space is 

zw„ -ek - Sii(k, iw„) -I]i2(k,iw„) \ / G'ii(k, iw„) Gi2(k, icj„) \ ^ / 1 

-i;2i(k,iw„) zcj„ + ek - S22(k,ia;„) y I G'2i(k,ia;„) G22{K'i^n) 1 \0 1 



(A2) 



where iw„ = 7r(2n+l)//3 are the fermionic Matsubara frequencies and /? = l/{k-BT) the inverse temperature. FormaUy 
the components of the Nambu-Green function can then be expressed as 

iu)„ + ek — S22(k, iujn) 



Gii(k,iaj„) = -G22(k, -ioJn) 



{iuJn - eic — Sll(k,ia)n))(i<.^n + Ek - 222(k, ioJn)) - Sl2(k, jtiJ„)S2l(k, itJn) 



^ /I • N ^ /I ■ \ Sl2(k, ioJn) 

Gi2(k,»t<;„) = G2i(k, — icjn) = ^- , (A3) 

(iuJn — £k - Sii(k,i<.i;„))(ja;n + £k — S22(k,iaj„)) - Si2(k,iaj„)S2i(k, icjn) 

where the d-dimensional dispersion is given by Sk = — 2t ^^ coa{kaaa)—n. The problem of calculatirig the components 
of the one-particle Nambu-Green function is then translated to the evaluation of the self-energies ||l^ . 

The equation of motion for the Nambu-Green function contains four-point correlation functions that cannot be 
evaluated analytically. In order to derive a perturbation expansion we add auxiliary source fields to the Hamiltonian 

and ultimately take the limit of these fields going to zero. This particular way of generating four-point correlation 
functions leads in the following to expressions for particle-particle scattering, which in turn leads to superconductivity 
in the Cooper channel. It is important to note that this approach does not generate particle-hole diagrams. 

The four-point correlation functions can now be written in terms of derivatives with respect to the auxiliary fields, 
and the equation of motion can be written 

/ -T-+tu'+l^ A(x)S(x-x') + U^4-^~ p*(x)\ , /l 0\ 

where A(x) = UGi2{x^ ,x) and A*(a;) = UG2i(x'^ ^ x). Here and in the following, we use the shorthand notation 
X = {1,t) for a space-time coordinate. Then, from Eq. (A5), the self-energy is given by E (a;, a;') = E"(x,x')-)-E^(a;, a;'), 
where 



and 







i:\x,a)g{a,x') = u[ ^ ^"^^ ]G{x,x') , (A7) 

V Sp'ix) " / 

where we used the notation f{a)g{a) = J daf{a)g{a). In the following, we also use the usual identity 6g{x,x') — 
g{x,a)S^{a,b)g^,x') §. 
^From Eq. ( |A7| ) we find 

E\x,x')^uJ2Msg{x,a)^-^^ , (A8) 

where the matrices Mi(2) are defined as 



Ml 



1 




M-2 





1 



(A9) 



and pi{x) = p{x),p2{x) ^ p*{x). 
We now introduce a vertex function 



r^^\x,x'\v)^u 



5Yi{x, x') 



(AlO) 



which obeys the equation 



off 



r'^-''\x,x'\y) = -U5{x-x')5{x^y)Mj 

+ US{x-x') \c{x,a)r^'\a,b\y)gib,x') 

+ J2Ms'g{x,a)r^-'Ha,b\y)g{b,c)r^''\c,x'\x) 

s' 

5r^''\a,x'\x) 



uY^Qix,} 



6ps{y) 



(All) 



where [...]off denotes the 2x2 matrix with zero diagonal elements and off-diagonal elements given by the matrix 
enclosed in brackets, and M""" is the transpose of the matrix M. The self-energy is given by 



E(a;, a;) = YP{x, x') + ^ g{x,a)r^'\a, x'\x) 



(A12) 



As Eq. ( All ) is not practically manageable, we wish to approximate it to a simpler form taking only ladder diagrams 
into account, which describes repeated scattering between two particles. This cont ribut ion is expected to dominate 
for small densities. A graphical analysis will easily show that the latter term in Eq. ( All ) corresponds to higher-order 
diagrams, so we omit it. The eight matrix elements of the vertex function then read 

rii{x,x'\y) = G2a{x,a)n^^^{a,b\y)Gp^f(b,c)r!:^^^{c,x'\x) , 

r^i^{x,x'\y) = US{x - x')Gio.ix,a)r^'^{a,b\y)Gf32{b,x') + G2aix,a)r^'^(a,b\y)Gp^(b,c)r^l\c,x'\x) , 

r^l\x,x'\y) = -U6{x - x')S{x - y) + U5{x - x')G2o.{x,a)r^^^{a,b\y)Gpi(b,x') 

+Gia{x,a)r^p{a,b\y)Gp^(b,c)r^\ {c,x'\x) , 
r^^{x,x'\y) = GiJx,a)ri'^{a,b\y)Gp^{b,c)r!^l\c,x'\x) , 
4f{x,x'\y) = -U6{x - x')S{x - y) + U5{x - x')Gi,(x, a)r^J(a, &|y)G^2(6, x') 

+G2a {x,a)n^il (a, b\y)G0^(b, c)r\^^ (c, x'\x) , 
4f {x,x'\y) ^U5{x ~ x')G2c.{x,a)r^^^p{aMv)Gpi(b,x') + Gi^{x,a)r'^^^p^ , 

where summation on repeated indices is understood. Similarly, the correction to the self-energy, E^, is given by 

i:.\^{x,x') = G2i{x,a)r{'J{a,x'\x) + G22{x,a)r!^]} {a,x'\x) , 
Ei„(x,x') = Gii{x,a)r'^^{a,x'\x) + Gi2{x,a)r!,^^{a,x'\x) . 

To further simplify this set of equations we expand it to second order in the anomalous Green's functions, Gi2(G2i), 
which yields the following approximate form of the vertex functions 



1 1-^ \x, X 

1 -^2 \X, X 

1 2]^ (X, X 

1 22 (X, X 

1 12 [X, X 

1 21 (a;, X 



''\y) « G22{x,a)4t\a,b\y)Gi2{b,c)r^l\c,x'\x) , 

■'\y) « USix ~ x')Gi2ix,a)4^^\a,b\y)Gi2(b,x') , 

%) « -US{x - x')6{x ~y) + US{x - x' )G 22 {x, a) r^l\a,b\y)Gn{b, x') 

''\y) « Gi2{x,a)4t\a,b\y)Gii{b,c)42\c,x'\x) , 

''\y) « -US{x ~ x')S{x -y) + US{x - x')Gn{x,a)4l\a,b\y)G22{b, x') 

■'\y) « US{x - x')G2i{x,a)r[^{a,b\y)G2i{b,x') , 



(A13) 



From Eqs. (A13) we conclude that r2i {x,x'\y) and rl^2 {^,x'\y) are related to the T-matrix equations above Tc 



m 

r2iix,x'\y) ^ r^l\x,x'\y) = -S{x ~ x')T{x,y) 
where T{x, y) in reciprocal space has the form 

U 



T{q, ie,n) = 



1 - Ux{q, ie.n) 



with 



Xili^Sr. 



jjs X! G22{k ~ q,iuj„ - iem)Gii{k,iuj„) 



N(i 



(A14) 



(A15) 



(A16) 



and lEn = Ittti / j3 a re the bosonic Matsubara frequencies. 

Using Eqs. (A14) and ( Al5| ), we obtain the self-energies to second order in 012(021): 

^\^{x,x') ^ G22{x,x')T{x' ,x) + G2i{x,a)G22{a,b)T(h,x)Gi2(b,x')T{x\x) , 

T,l2{x,x') = G22ix,a)Gi2ia,b)T(b,x)Gii{b,x')T{x',a) , 

^liix.x') = Gn{x,a)G2i{a,b)T{b,x)G22{b,x')T(x',a) , 

1122(2^; 2;') ~ Gii{x,x')T{x' ,x) + Gi2{x,a)Gii{a,b)T(b,x)G2i(b,x')T{x',x) . 



(A17) 



This expansion is v alid forA/W^C 1, where W = 2dt is the bandwidth. 

To solve Eqs. (A3) and ( A15 )-( A17 ) one would also have to fix the chemical potential from the particle number 
using 



P{T, m) = lin\ -3T7 y2 '^(^' *^") exp(zw„?7) , 



(A18) 



where p is t he electron concentration per spin and is defined in the interval [0,1]. Thus, the set of Eqs. ( |A3| ) and 
( A15| )-(A18) represents a set of nonlinear self-consistent equations, which must be solved numerically. 

We note that an expansion of the final equations, Eqs. ( A13| )-(A17), to first order in U simply yields the well- 
known BCS expressions. To second order in [/, the result is identical to that of Martin-Rodero and Flores |l3]. The 
second-order expansion was found to yield the same gap equation as in BCS, but with a renormalized interaction. 



APPENDIX B: DERIVATION OF MOMENT EQUATIONS 

Introduce the diagonal and ofF-diagonal one-particle Green's function 

and 

where T,- is the time-ordering operator, as well as the associated spectral functions 

A{k, to) = IniG(fc, to + iS) , 

n 

B{k, oj) = IniF(fc, uj + iS) . 

n 

To construct approximate expressions for the spectral functions, it is useful to consider the frequency moments 

An{k) 

Bn{k) 
With the d-dimensional dispersion defined as 

e-k - 

the exact four first diagonal moments are given by ||9[] 

Ao{k) = 1 

Ai{k) ^ ai = eu - ^JL - pU 
A2{k) =02 = (ffc - ^J-f - 2(efc - ^Ji)pU + pU^ 

Asik) = 03 = (efc - t^f - 3(efc - ti)pU + (e^ - p)U^p{2 + p) + (e^ - ti)p''U^ 
-U^p+U^p{l-p)K{k) , 



dujuj''B{k,uj) . 

d 

-2t'^cos{kin) , 



where 



p{l~p)K{k)=K + KUk) 



KUk) = ^Y.^^Je'''^'^^-'^'H{n^<,n,^) ~ p" 



„t J 



\pjaC,gCifjCi„) {Cj^CiffCja-Cia-} f 



and the first two off-diagonal moments are given by 

Bo{k) = 



(Bl) 
(B2) 

(B3) 
(B4) 

(B5) 
(B6) 

(B7) 



(B8) 

(B9) 

(BIO) 

(Bll) 

(B12) 



(B13) 
(B14) 

(B15) 



B,{k) = -UJ2{c^lC^^) 



(B16) 
(B17) 



In the normal state the first four frequency moments have recently been used in conjunction with an ansatz of 
the form A{k,u;) = ai{k)S{iLj — ^i(fc )) -|- a2{k)S{ui — ^2{k))- The fc-dependence in Eq. (B13) is approximated by its 
momentum average. Only Eq. (B14) contributes to the average and is given by [|| 



p{l - P)^ = -^ E E ".(fc)efc"/(a(fc)) 

(B18) 



N 

k i 



where nf{uj) is the Fermi distribution function. The chemical potential is obtained by fixing the density using 

k,i 

To analyze the superconducting state we make the assumption that the general structure of the diagonal part of the 
self-energy does not change in the superconducting state. As the diagonal part of the self-energy is related primarily 
to pair-interaction physics, we expect this to be a reasonable approximation. Thus, we express the self-energy as 

1 ^ "iW , "2(fc) ,„ . 

z-efe + S|(fc,z) z ~ ni{k) z - n2ik) ' ^ ' 

where ai(fc), a2{k), Qi(k), Q2{k) are functions to be determined. Note that this approximation neglects lifetime 
effects. 

The diagonal and off-diagonal Green's function in the superconducting state is now given by Dyson's equation 

r,(k A = z + ek-^i{k,z) 

^^'' {z-ek + ^^{k,z)){z + ek-^i{k,z))^\Ak\^ ^''^'> 

^^^'^^" (z-efe + ST(fc,z))(z + efc + I]x(fc,z))-|Afe|2 ' ^^^^^ 

where Sj^(/c,z) — — E|(fc,— z). In Dyson's equation we have made the additional assumption that the off-diagonal 
part of the self-energy, A^, is independent of frequency. This assumption is likely to be wrong for large coupling 
strengths, where A^ might contain additional structure related to the existence of local pairs. 
Combining the above, we can write Green's functions for frequencies on the real axis as 



ai{k){Lo + ni{k)){ Lo^ ~ nljk)) + a2{k){uj + n 2{k)){uo^ - nl{k)) 

{^^~u:l{k)){u:^~ul{k)) 
Afc[Lj ^ ~ (ai(fc)r>2(fc) + a2{k)ni {k)f] 
{u;^-u;Uk)){u;^-u;Uk)) 

where ±wo(fc) and ±wi(fc) are the dispersions of the resulting four poles given by 

'^0,i('=) = ^ (^fe + ^i^'') + ^2(*=) ^ Vl^fc + ^?(^) + ^U^^)? - 4p?(A;)n2(fc) + Al{ai{k)n2ik) + a2(k)ni{k))2]\ . (B25) 



Gt (k, -) = ^^-fM(^J^^^^fP^ ^ ^ 

F{k,cu)^ ^''^ ~,r'''^^n't\'^7^'i^n^'^''^^'^ , (B24) 



The corresponding spectral functions, Eqs. (B3) and (B4), can then be written as A(k,Lu) = ai{k)S{Lu — loq) + 



a2ik)S{Lu + LUo) +a3{k)S{uj ~uji) + a4{k)S {lu + lui) and B{k,uj) — (3i(5(uj — ujq) — 6{lo + loo)) + I32{5{oj — loi) —5{uj + uji)). 
The weights are found as the residues of the poles and given by 

- ,,, ai(fc)(c^o(fc) + »i(fc))K(fc)^ - n2{kf) + «2(fc)(^o(fc) + n2{k)){Lo^{kf - n^jkf) 

"^^'^ = 2u:,mMk?-^.{k?) ^""''^ 

. ai(fc)(-a;o(fc) + n^{k)){u:o{kf - n2{kf) + a2(fc)(-c.o(fc) + n2{k)){uJo{kf - n,{kf) 

"^^'^ -2uo{k){Mkr ~ u;,{kr) ^^''^ 

- ,,, ai(fc)(a;i(fc) + »i(fe))(u;i(fc)^ - ^^W") + «2(fc)(^i(fc) + 1^2(fc))K(fc)' - ^i(fc)') ,„„„, 

"^^'^ ^ 2umi-,ikr-Mkr) ^""''^ 



ai(fc)(-a;i(fc) + r>i(fc))(a>i(fc )2 - 1^2^) + a2(fc)(-c^i(fc ) + r!2(fc))K(fc)2 - r!i(fc)2) 

Akitooik)^ ~ {Mk)^2{k) + a2 {mi{k)f' 

2wo(fc)(wo ^^l) 
Ak{uJo{k? - (ai(fc)1^2(fc) + a2(fc)r!i(fc))2^ 



«4(fc) = ....,.....)2_,^(fc)2) (B29) 
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Knowing the spectral functions, we now insert them into the moment equations (B8)-(B12) and obtain the foUowing 
for the superconducting state: 



ai{k) + a2{k) — 1 
ai{k)ni{k) + a2{k)n2ik) = ai 
ai{k)nl{k) + a2{k)nl{k) = as - A^ 
ai{k)nl{k) + a2{k)nl{k) = 03 - A^oi 



(B32) 
(B33) 
(B34) 
(B35) 



The first two off-diagonal moment equations, Eqs. ( p316| ) and (B17), are automatically satisfied by our Dyson's 
equation. 

The resulting moment equations for the superconducting state can be solved and uniquely determine the form of 
the spectral functions. The second off-diagonal moment equation yields the gap equation 



U 



/oc 
duj 7if{Lj)B{k^ uj) 
-00 



(B36) 



We immediately see that in this approach the order parameter is a constant, A^ = A. We now have to evaluate 
t he or der p ar amet er , chem ical potential, band correction, ai{k) , a2{k) , fli{k) , and r22(fc) self-consistently, using Eqs. 
(p8|)-(pT2|), (pl8|), (|bT9|), and (|B36|). 
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FIG. 1. Diagonal one-particle spectral function, A(k,a;) vs. u) for various momenta along the diagonal of the Brillouin zone, 
k = (n, n)7r/16, for U/t = 6.0, T/t = 0.05, and p = 0.1. We have used an external damping of 5/t — 0.1. After self-consistent 
calculation of the coupled nonlinear equations, we obtain the chemical potential, fi/t — —4.438, and the order parameter, 
A/t = 1.04. 
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FIG. 2. Off-diagonal one-particle spectral function, i3(k,a;) vs u for various momenta along the diagonal of the Brillouin 
zone. Same parameters as in Fig. hi. 
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FIG. 3. Real part of the self-energy, Re[E(k, ui)] vs to for various momenta along the diagonal of the Brillouin zone. Same 
parameters as in Fig. nl. 
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FIG. 4. — Im[G2(q,tj)] vs oj for different momenta along the diagonal of the Brillouin zone, q — {m,m)n/16. Same 
parameters as in Fig. 
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FIG. 5. Quasiparticle dispersion for the two poles in the normal state at T/t = 0.2 taken along the diagonal in the Brillouin 
zone for U/t = 4, p = 0.1. 
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FIG. 6. Amplitudes of the two poles in the normal state at T /t — 0.2 taken along the diagonal in the Brillouin zone for 
U/t = 4, p = 0.1. 
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FIG. 7. Quasiparticle dispersion for the four poles in the superconducting state at T/t = 0.05 taken along the diagonal in 
the Brillouin zone for U/t = 4, p = 0.1. For comparison, T-matrix results are included as symbols. 
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FIG. 8. Amplitudes for the four poles in the superconducting state at T jt — 0.05 taken along the diagonal in the Brillouin 
zone for U/t — 4, p = 0.1. For comparison, T- matrix results are included as symbols. 
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